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PREFACE 


The  study  reported  herein  was  authorized  by  the  Office,  Chief  of  Engi¬ 
neers  (OCE) ,  US  Army  Corps  of  Engineers,  and  funded  through  the  Coastal  Engi¬ 
neering  Functional  Area  of  Civil  Works  Research  and  Development,  under  Work 
Unit  31616.  The  project  was  monitored  by  Messrs.  John  H.  Lockhart,  Jr.,  and 
John  G.  Housley,  OCE  Technical  Monitors.  The  study  was  conducted  at  the 
Coastal  Engineering  Research  Center  (CERC)  of  the  US  Army  Engineer  Waterways 
Experiment  Station  (WES).  Dr.  C.  Linwood  Vincent,  CERC,  is  Program  Manager  of 
the  Coastal  Engineering  Functional  Area. 

This  report  was  prepared  by  Mr.  John  P.  Ahrens,  Research  Oceanographer, 
Wave  Research  Branch  (CW-R) ,  Wave  Dynamics  Division  (CW),  CERC.  Assisting 
Mr.  Ahrens  in  conducting  the  study  were  the  following  CERC  employees: 

Ms.  Karen  P.  Zirkel  and  Messrs.  Louis  Myerele  and  Martin  F.  Titus,  Engineering 
Technicians;  Messrs.  John  Heggins,  Computer  Assistant,  and  Leland  Hennington, 
Summer  Aide,  who  helped  to  analyze  the  data;  and  Eng.  Gisli  Viggosson  on  temp¬ 
orary  assignment  from  the  Icelandic  Harbour  Authority,  Reykjavik,  Iceland. 

Work  was  performed  under  direct  supervision  of  Messrs.  D.  D.  Davidson, 
CW,  and  C.  Eugene  Chatham,  Chief,  CW;  and  under  general  supervision  of 
Dr.  James  R.  Houston  and  Mr.  Charles  C.  Calhoun,  Jr.,  Chief  and  Assistant 
Chief,  CERC,  respectively.  This  report  was  edited  by  Ms.  Shirley  A.  J. 
Hanshaw,  Information  Products  Division,  Information  Technology  Laboratory, 

WES. 

Commander  and  Director  of  WES  during  publication  of  this  report  was 
COL  Dwayne  G.  Lee,  CE.  Technical  Director  was  Dr.  Robert  W.  Whalin. 
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CHARACTERISTICS  OF  REEF  BREAKWATERS 


PART  I :  INTRODUCTION 

1.  A  reef  breakwater  is  a  low-crested  rubble-mound  breakwater  without 
the  traditional  multilayer  cross  section.  This  type  of  breakwater  is  little 
more  than  a  homogeneous  pile  of  stones  with  individual  stone  weights  similar 
to  those  ordinarily  used  in  the  armor  and  first  underlayer  of  conventional 
breakwaters. 

2.  In  recent  years  a  number  of  low-crested  breakwaters  have  been  built 
or  considered  for  use  at  a  variety  of  locations.  Most  of  these  structures  are 
intended  to  protect  a  beach  or  reduce  the  cost  of  beach  maintenance.  Other 
applications  include  protecting  water  intakes  for  power  plants  and  entrance 
channels  for  small-boat  harbors  and  providing  an  alternative  to  revetment  for 
stabilizing  an  eroding  shoreline.  In  situations  where  only  partial  attenua¬ 
tion  of  waves  on  the  leeside  of  a  structure  is  required,  or  possibly  even 
advantageous,  a  low-crested  rubble-mound  breakwater  is  a  logical  selection. 
Since  the  cost  of  a  rubble-mound  breakwater  increases  rapidly  with  the  height 
of  the  crest,  the  economic  advantage  of  a  low-crested  structure  over  a  tradi¬ 
tional  breakwater  that  is  infrequently  overtopped  is  obvious.  Because  the 
reef  breakwater  represents  the  ultimate  in  design  simplicity,  it  could  be  the 
optimum  structure  for  many  situations.  Unfortunately,  the  performance  of  low- 
crested  rubble-mound  structures,  particularly  reef  breakwaters,  is  not  well 
documented  or  understood. 


Background 

3.  A  number  of  papers  have  noted  that  armor  on  the  landside  slope  of  a 
low-crested  breakwater  is  more  likely  to  be  displayed  by  heavy  overtopping 
than  armor  on  the  seaward  face  (Lording  and  Scott  1971,  Raichlen  1972,  and 
Lillevang  1977).  Raichlen  discusses  the  characteristics  of  overtopping  over 
the  c rest  and  the  inherent  complexity  of  the  problem.  Walker,  Palmer,  and 
Dunham  (1975)  give  a  carefully  reasoned  discussion  of  the  many  factors  influ¬ 
encing  stability  of  heavily  overtopped  rubble-mound  breakwaters.  They  also 
show  a  figure  which  suggests  what  armor  weight  is  required  for  stability 


on  the  backside  of  a  low-crested  breakwater.  Unfortunately,  the  data  scatter 
shown  in  the  figure  undermines  confidence  in  the  suggested  armor  weights. 

4.  In  Australia,  the  breakwater  at  Rosslyn  Bay  was  damaged  severely 
during  Cyclone  David  in  1976  (Bremner  et  al.  1980).  The  crest  height  of  the 
structure  was  reduced  as  much  as  4  m  but  still  functioned  effectively  as  a 
submerged  breakwater  for  over  2  years  until  it  was  repaired.  Based  on  the 
surprisingly  good  performance  of  the  damaged  Rosslyn  Bay  breakwater  and  the 
findings  from  model  tests,  a  low-crested  design  was  chosen  for  the  breakwater 
at  Townsville  Harbor,  Australia.  This  breakwater  is  unusual  because  it  was 
built  entirely  of  stone  in  the  3-  to  5-ton*  range  (Bremner  et  al.  1980).  Reef 
breakwaters,  as  described  in  this  paper,  are  very  similar  to  the  Townsville 
breakwater  except  a  wider  gradation  of  stone  was  used  in  the  model  breakwater 
tests  discussed  herein. 

5.  Seelig  (1979)  conducted  an  extensive  series  of  model  tests  to  deter¬ 
mine  wave  transmission  and  reflection  characteristics  of  low-crested  break¬ 
waters,  including  submerged  structures.  From  these  tests  Seelig  concluded 
that  the  component  of  transmission  resulting  from  wave  overtopping  was  very 
strongly  dependent  on  the  relative  freeboard  (i.e.,  freeboard  divided  by  inci¬ 
dent  significant  wave  height).  Recent  work  by  Allsop  (1983)  with  multi¬ 
layered,  low-crested  breakwaters  shows  that  wave  transmission  is  strongly 
dependent  on  a  dimensionless  freeboard  parameter  which  includes  the  zero¬ 
crossing  period  of  irregular  wave  conditions.  Allsop  did  not  find  substantial 
wave  period  dependency  in  his  evaluation  of  breakwater  stability.  He  indi¬ 
cates,  however,  that  since  wave  transmission  (which  largely  results  from  over¬ 
topping)  is  dependent  on  period,  then  possible  stability  of  the  backside  slope 
would  also  be  a  function  of  wave  period. 

Scope 

6.  A  study  currently  being  conducted  at  the  US  Army  Engineer  Waterways 
Experiment  Station's  Coastal  Engineering  Research  Center  is  intended  to  docu¬ 
ment  the  performance  of  low-crested  breakwaters.  This  paper  discusses  labora¬ 
tory  model  tests  of  reef  breakwaters  and  provides  information  on  their  stabil¬ 
ity  to  wave  attack,  wave  transmission  and  reflection  characteristics,  and  wave 
energy  dissipation. 


PART  II:  LABORATORY  SETUP  ANT)  TECHNIQUES  USED 


7.  To  date,  205  two-dimensional  laboratory  tests  of  reef  breakwaters 
have  been  completed.  These  tests  were  conducted  in  a  61-cm-wide  channel 
within  CERC's  1.2-  by  4.6-  by  42.7-m  tank  (Figure  1).  All  tests  were 
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Figure  1.  Plan  view  of  wave  tank  and  test  setup 

conducted  with  irregular  waves.  The  spectra  used  had  wave  periods  of  peak 
energy  density  T^*  ranging  from  about  1.45  to  3.60  sec,  and  water  depth  at 
the  structure  dg  ranged  from  25  to  30  cm.  Signals  to  control  the  wave  blade 
were  stored  on  magnetic  tape  and  transferred  to  the  wave  generator  through  a 
computer  data  acquisition  system  (DAS).  For  this  study  four  files  were  stored 
on  the  tape  which  could  produce  a  spectrum  with  a  distinct  period  of  peak 
energy  density.  Table  1  gives  the  nominal  period  of  peak  energy  density  for 
each  file. 

8.  If  there  were  no  attenuation  of  the  signal  to  the  wave  generator, 
the  files  used  were  intended  to  produce  a  saturated  spectrum  at  all  frequences 
above  the  frequency  of  peak  energy  density  for  the  water  depth  at  the  wave 
blade.  For  frequencies  lower  than  those  of  the  peak,  the  energy  density  de¬ 
creased  rapidly.  This  procedure  produced  a  spectrum  of  the  Kitaigorodskii 
type  as  described  by  Vincent  (1981).  The  amplitude  of  the  signal  to  the  wave 
generator  was  attenuated  by  a  10-turn  potentiometer  in  a  voltage  divider 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  C) . 
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Table  1 


H 


Period  of  Peak  Energy  Density 


for  Each  Tape  File 


Approximate 
T  ,  sec 


network  which  allowed  control  of  the  wave  heights  generated.  In  addition,  the 
waves  were  generated  in  a  water  depth  25  cm  greater  than  at  the  breakwater  and 
shoaled  to  the  water  depth  at  the  structure  over  a  1-V  on  15-H  slope  (see 
Figure  2) .  This  setup  ensures  that  severe  conditions  can  be  developed  at  the 
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Figure  2.  Cross  section  of  test  channel 

structure  site.  Incident  zero-moment  wave  heights  ranged  from  about  1 

to  18  cm. 

9.  Three  parallel  wire-resistance  wave  gages  were  used  in  front  of  the 
breakwater  to  resolve  the  incident  and  reflected  wave  spectra  using  the  method 
of  Goda  and  Suzuki  (1976),  and  two  wave  gages  were  placed  behind  the  structure 
to  measure  the  transmitted  wave  height.  The  location  of  gages  is  shown  in 
Figure  2.  During  data  collection  gages  were  sampled  at  a  rate  of  16  times  per 
second  for  256  sec  by  the  same  DAS  which  controlled  the  wave  generator  motion. 
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10.  Two  types  of  model  tests  were  conducted  during  this  study:  stabil¬ 
ity  and  previous  damage  tests.  Each  type  followed  a  prescribed  sequence. 

Stability  Tests 

11.  For  a  stability  test  the  following  test  sequence  was  used: 

a.  Rebuild  the  breakwater  from  the  previously  damaged  condition. 

b.  Survey  the  breakwater  to  document  its  initial  condition. 

c.  Calibrate  the  wave  gages. 

d.  Select  the  tape  file  and  signal  attenuation  setting. 

e.  Start  the  wave  generator  and  run  waves. 

Collect  wave  data  (several  or  more  times). 

£.  Stop  the  wave  generator. 

h.  Survey  the  breakwater  to  document  its  final  condition. 

The  duration  of  wave  action  was  from  1-1/2  hr  for  a  test  using  the  File  1 
spectrum  to  3-1/2  hr  for  a  File  4  spectrum.  Generally,  the  technicians 
observing  the  tests  thought  that  most  of  the  stone  movement  occurred  during 
the  first  10  or  15  min  of  wave  generation,  so  the  final  survey  is  regarded  as 
an  equilibrium  profile  for  the  structure.  In  rebuilding  the  breakwater  the 
technicians  rarely  touched  the  stone  but  merely  pushed  it  around  by  foot  until 
the  shape  conformed  to  the  desired  initial  profile.  This  procedure  was  a  con¬ 
scious  effort  to  avoid  overly  careful  placement  of  the  stone.  Outlines  of  the 
desired  initial  profile  were  fixed  to  the  walls  of  the  testing  channel,  and  a 
moveable  template  was  used  to  ensure  that  the  initial  profile  was  reasonably 
close  to  the  desired  profile.  Initial  configuration  of  the  breakwater  for  a 
stability  test  was  a  narrow,  trapezoidal  shape  with  seaward  and  landward 
slopes  of  IV  on  1.5H  (Figure  3).  Crest  widths  were  three  typical  stone  dimen¬ 
sions  wide,  using  the  cube  root  of  the  volume  of  the  median  weight  stone 
as  the  typical  dimension  d^^  .  Figure  3  also  shows  a  typical  profile  after 
moderately  severe  wave  attack  during  a  stability  test.  Wave  transmission  and 
reflection  also  were  measured  during  a  stability  test. 
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Figure  3.  Cross-sectional  view  of  initial  and  typical  damaged  reef  profiles 

(swl  denotes  still-water  level) 

Previous  Damage  Tests 

12.  Previous  damage  tests  were  conducted  to  answer  the  question  of  how 
the  breakwater  would  perform  for  moderate  wave  conditions  after  it  had  been 
damaged  by  very  severe  wave  conditions.  For  previous  damage  tests  there  was 
very  little  readjustment  of  the  damage  profile  from  test  to  test;  conse¬ 
quently,  the  breakwater  was  not  rebuilt  at  the  end  of  a  test.  No  stability 
information  was  obtained  from  these  tests,  and  the  duration  of  wave  action  was 
only  half  an  hour;  however,  wave  transmission  and  reflection  were  measured. 
Previous  damage  tests  were  performed  in  the  following  sequence: 

a.  Survey  breakwater  for  last  test  which  becomes  initial  survey 
for  current  test. 

b.  Calibrate  wave  gages. 

c.  Select  wave  file  and  signal  attenuation  setting. 

d.  Start  generator  and  run  waves  for  half  an  hour. 

e.  Collect  wave  data  (two  or  three  times). 

£.  Stop  wave  generator. 

£.  Survey  breakwater  as  noted  above  in  Step  1. 

13.  All  205  of  the  completed  tests  of  this  study  can  be  divided  logi¬ 
cally  into  10  subsets  or  test  series.  Because  of  the  test  plan,  stability 
test  series  have  odd  numbers,  and  previous  damage  test  series  have  even  num¬ 
bers.  Table  2  lists  the  basic  information  about  each  subset. 

14.  Two  different  sizes  of  stone  were  used  during  this  study.  For  sub¬ 
sets  I  through  6  an  angular  quartzite  with  a  median  weight  of  17  g  was  used. 


Table  2 


Basic  Data  for  Each  Subset 


Subset 

No. 

No . 
of 

Tests 

Water 

Depth 

d  ,  cm 
s 

G  rest 

Height 

"as  built" 

h '  ,  cm 
c 

Median 

Stone  Weight 

W50  ’  g 

Area  of 

Breakwater 

Cross  Section 

.  2 

At  ,  cm 

1 

27 

25 

25 

17 

1,170 

2 

3 

25 

NA* 

17 

1,170 

3 

29 

25 

30 

17 

1,560 

4 

12 

25 

NA 

17 

1,560 

5 

41 

25 

35 

17 

2,190 

6 

1 1 

25 

NA 

17 

2,190 

7 

38 

25 

32 

71 

1 ,900 

8 

26 

25 

NA 

71 

1,900 

9 

13 

30 

32 

71 

1,900 

10 

5 

30 

NA 

71 

1 ,900 

*  NA  denotes  not  applicable  to 

previous  damage 

test  series. 

and  for  subsets  7  through  10  a  blocky  to  angular  diorite  with  a  median  weight 
of  71  g  was  used.  Photo  1  depicts  the  stone,  and  Table  3  summarizes  informa¬ 
tion  about  it. 


Table  3 

Stone  and  Gradation  Characteristics 


Characteristic 

Quartzite 

Diorite 

27  weight  (g) 

7.0 

14.0 

Median  weight,  W  ^  (g) 

17.0 

71.0 

98%  weight  (g) 

28.0 

139.0 

2 

Density  (g/cm  ) 

2.63 

2.83 

Porosity  (7) 

45 

44 

Profile  Surveys 

15.  Initial  and  final  profiles  of  the  reef  were  obtained  by  survey. 

The  survey  rods  had  feet  attached  with  ball-and-socket  connectors.  For  the 
small  stone  used  for  subsets  1  through  6,  the  foot  of  the  survey  rod  had  a 
diameter  of  2.54  cm;  and  for  the  somewhat  larger  stone  used  in  subsets  7 
through  10,  the  foot  of  the  survey  rod  had  a  diameter  of  3.81  cm.  Three  pro¬ 
files  were  used  to  establish  an  average  profile  for  the  reef.  One  profile 
line  was  exactly  in  the  center  of  the  wave  channel,  and  the  other  two  profile 
lines  were  15  cm  on  either  side  of  center.  The  survey  interval  along  the 
channel  was  3.05  cm. 


PART  III:  STABILITY  AND  PERFORMANCE  RESULTS 


16.  The  report  herein  consolidates  findings  from  all  of  the  data  sub¬ 
sets  identified  in  Table  2  into  general  conclusions  about  the  stability  and 
performance  characteristics  of  reef  breakwaters.  Specific  characteristics  in¬ 
clude  the  stability  of  reef  breakwaters  to  irregular  wave  attack,  wave  trans¬ 
mission  over  and  through  the  breakwater,  wave  reflection  from  the  breakwater, 
and  dissipation  of  wave  energy.  A  mathematical  model  is  developed  for  each 
characteristic  which  provides  a  simple  method  to  summarize  findings  from  this 
study  and  a  convenient  way  to  furnish  results  to  potential  users.  These  math¬ 
ematical  models  are  intended  to  work  together  with  the  stability  model  fur¬ 
nishing  the  equilibrium  crest  height  to  both  transmission  and  reflection 
models  which  together  are  used  to  estimate  the  amount  of  energy  dissipated  by 
the  reef. 

Stability  to  Irregular  Wave  Attack 

17.  The  stability  of  reef  breakwaters  will  be  quantified  by  damage  or 
lack  of  damage  during  a  test,  the  most  important  aspect  of  which  is  the  reduc¬ 
tion  in  crest  height  caused  by  wave  attack.  This  aspect  of  stability  is 
important  because  the  performance  of  a  reef  breakwater  will  be  judged  largely 
on  its  wave  transmission  characteristics.  Wave  transmission  is  very  sensitive 
to  crest  height  relative  to  water  level. 

Crest  height  reduction  factor 

18.  One  of  the  most  effective  methods  to  evaluate  damage  to  a  reef 

breakwater  is  to  use  the  ratio  of  the  crest  height  at  the  completion  of  a  test 

to  the  height  at  the  beginning  of  the  test  before  waves  have  been  run.  This 

ratio,  h  /h'  ,  will  be  referred  to  as  the  crest  height  reduction  factor.  For 
c  c 

comparing  damage  within  a  subset,  h^/h^  is  effective  because  it  inherently 
accounts  for  the  random  variation  of  one  to  two  centimeters  in  the  constructed 
crest  height  from  test  to  test  within  a  subset.  Another  advantage  of  the 
crest  height  reduction  factor  is  that  all  stability  subsets  have  the  same 
natural  limiting  values  of  1.0  and  0.0. 

Stability  number  and 

spectral  stability  number  comparison 


variables  to  explain  damage  would  be  one  similar  to  the  stability  number  used 
by  Hudson  and  Davidson  (1975).  The  following  definition  is  used  for  the 
stability  number  for  tests  with  irregular  waves: 


where  w  is  the  density  of  stone  and  w  is  the  density  of  water.  Since 

r  W  3 

these  tests  were  conducted  in  fresh  water,  w  *1.0  g/cm  .  As  far  as  the 

w 

stability  tests  of  reef  breakwaters  are  concerned,  it  was  apparent  that  tests 
with  a  higher  period  of  peak  energy  density  did  more  damage  than  similar  tests 
with  a  smaller  period  of  peak  energy  density.  This  finding  is  consistent  with 
the  results  of  a  study  conducted  by  Gravesen,  Jensen,  and  Sorensen  (1980)  on 
the  stability  of  high-crested,  rubble-mound  breakwaters  exposed  to  irregular 
wave  attack.  According  to  the  stability  analysis  of  Gravesen,  the  spectral 
stability  number  is  defined 


where  L  is  the  Airy  wave  length  calculated  using  T  and  the  water  depth 
P  P 

at  the  toe  of  the  reef  d 

s 

20.  Figures  A  through  8  show  comparisons  of  the  effectiveness  of  the 
stability  number  and  the  spectral  stability  number  in  accounting  for  damage  to 
reef  breakwaters.  In  Figures  A,  5,  6,  7,  and  8  the  crest  height  reduction 
factor  is  plotted  versus  the  traditional  stability  number  and  the  spectral 
stability  number  for  stability  subsets  I,  3,  5,  7,  and  9,  respectively.  The 
figures  show  that  there  is  less  scatter  in  the  damage  trends  when  they  are 
plotted  versus  the  spectral  stability  number.  They  also  show  that  there  is 
little  or  no  damage  for  spectral  stability  numbers  less  than  about  six  but 
that  damage  increases  rapidly  for  spectral  stability  numbers  above  eight.  Tn 
the  following  analysis  the  spectral  stability  number  will  be  used  to  define 
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Crest  height 

Figure  8. 


reduction  factor  versus  the  spectral 
stability  number 

Stability  comparisons  for  subset  9 


the  relative  severity  of  wave  attack  on  reef  breakwaters. 
Secondary  stability  factors 


21.  Data  analysis  and  observation  of  the  laboratory  tests  indicate  that 
several  factors  other  than  the  spectral  stability  number  have  a  quantifiable 
influence  on  the  stability  of  reef  breakwaters.  Figure  9  will  help  identify 
what  will  be  referred  to  as  secondary  stability  factors  or  variables.  In 
Figure  9  the  damage  trends  for  all  five  stability  subsets  are  shown  using  sub¬ 
jectively  drawn  curves.  Figure  9  shows  the  relative  crest  height  h  /d  (see 

c  s 

Figure  3)  as  a  function  of  the  spectral  stability  number.  For  intercomparing 

damage  trends  between  subsets,  the  variable  h  /d  is  better  than  h  /h'  . 

c  s  c  c 

When  various  subsets  are  plotted  using  IWh^  »  the  data  trends  tend  to  fall 
on  top  of  each  other,  especially  for  N*  <  8  .  Using  h^/d^  to  show  damage 
trends  spreads  the  data  out  so  that  subsets  can  be  distinguished  and  provides 
better  orientation  by  showing  the  swl. 

22.  Relative  exposure  to  wave  action.  One  secondary  stability  factor 
is  the  relative  exposure  of  the  structure  to  wave  action.  Submerged  break¬ 
waters  are  much  less  exposed  to  wave  attack  than  breakwaters  with  crests  above 
the  water  level.  Water  overlying  a  submerged  crest  greatly  dampens  wave 
impact  forces  and  attenuates  the  lift  and  drag  forces  on  the  stone.  This  fac¬ 
tor  is  illustrated  in  Figure  9  where  structures  with  the  greater  initial  rela¬ 
tive  height  hVdg  have  their  height  reduced  more  rapidly  with  increasing  N* 
than  structures  with  lower  initial  relative  height.  In  Table  4,  which  can  be 
used  with  Figure  9  to  evaluate  the  influence  of  secondary  stability  factors, 
the  average  value  of  initial  relative  crest  height  hVdg  is  given  by  subset 
along  with  two  other  secondary  stability  factors,  the  bulk  number  and  the  "as 
built"  effective  reef  slope  C'  ,  which  are  discussed  below.  Subsets  I  and  5, 
which  represent  tests  using  the  same  stone  size  and  water  depth,  illustrate 
the  influence  of  l//dg  on  stability.  Figure  9  shows  that  the  wide  dif¬ 
ference  in  initial  relative  height  of  these  structures  is  maintained  until  N* 

s 

is  about  6.0;  however,  when  noticeable  stone  movement  starts  at  about  N* 

s 

=  6  ,  the  difference  in  relative  heights  for  the  reefs  of  the  two  subsets 
tends  to  decrease  with  increasing  value  of  N*  .  For  the  most  severe  condi¬ 
tions  at  about  N*  =  17  ,  the  difference  in  relative  height  between  the  two 
subsets  is  not  very  large.  Based  on  analysis  of  all  the  data,  it  is  concluded 
that  the  greater  the  initial  height  of  the  reef  the  more  vulnerable  it  is  to 
wave  attack. 
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23.  Influence  of  reef  bulk.  Subsets  1  and  5  can  be  used  also  to 
illustrate  the  influence  of  size  or  bulk  of  the  reef  on  stability.  Even 
though  the  difference  in  relative  height  for  the  two  subsets  narrows  with 
increasing  N*  ,  the  crest  heights  of  the  reefs  of  subset  5  always  are  higher 
than  those  of  subset  1.  In  fact.  Figure  9  shows  that  the  relative  position  of 
the  trends  for  subsets  I,  3,  and  5  are  maintained  such  that  the  larger  struc¬ 
ture  always  has  a  greater  crest  height  than  the  smaller  structure  for  a  given 
value  of  N*  .  In  order  to  intercompare  the  stability  of  all  subsets,  a 
general  measure  of  breakwater  size  is  needed  which  will  be  consistent  with  the 
data  trends  shown  in  Figure  9.  Within  this  context,  the  variable  which  best 
characterizes  the  size  of  the  reef  breakwater  is  called  the  bulk  number  B 

n 

and  is  defined  as 


W  \2 
— 

Wr/ 


where 

2 

A  «*  area  of  breakwater  cross  section,  cm 

c  3 

w^_  =  unit  weight  of  stone  g/cm 

d<.Q  =  dimension  of  stone,  cm 

24.  Bulk  number  can  be  described  as  the  equivalent  number  of  median 
stones  per  median  stone  width  in  the  breakwater  cross  section.  Equivalent  is 
used  because  B^  does  not  include  the  influence  of  porosity  which  is  about 
45  percent  for  the  two  stone  gradations  used  in  this  study.  The  value  of  the 
bulk  number  lies  in  its  ability  to  explain  the  rather  straightforward  behavior 
of  the  relative  location  of  the  damage  trends  for  subsets  1,  3,  and  5  in  Fig¬ 
ure  9.  It  also  explains  the  rather  anomalous  behavior,  such  as  that  of  the 
trend  for  subset  9  crossing  the  trend  for  subset  1.  At  first  it  seems  sur¬ 
prising  that  the  reefs  of  subset  9  degrade  faster  than  those  of  subset  1,  con¬ 
sidering  that  the  reefs  of  subset  9  have  the  greater  cross-sectional  area  (see 
Table  2).  However,  when  the  bulk  number  is  used  to  measure  the  size  of  the 
reef  rather  than  the  cross-sectional  area,  the  relative  behavior  of  the  damage 
trends  for  subsets  1  and  9  seems  more  plausible.  Subsets  1  and  9  have  bulk 
numbers  of  337  and  222,  respectively,  indicating  that  the  reefs  of  subset  1 


Uyiy 

3' 


EE? 


m 


K?' 


8 


•S'd.v.h'. 


>  -'YV* 


have  more  stone  in  the  cross  section  than  the  reefs  of  subset  9.  All  the  data 
appear  to  indicate  that  when  the  relative  severity  of  wave  attack  is  based  on 
the  spectral  stability  number  the  stability  of  the  reef  correlates  better  with 
the  number  of  stones  in  the  cross  section  than  with  the  absolute  size  of  the 
cross  section.  Other  factors  being  equal,  a  reef  with  a  large  bulk  number  is 
more  stable  than  a  reef  with  a  small  bulk  number  because  there  are  more  stones 
to  dissipate  wave  energy  and  to  shelter  other  stones  from  wave  forces. 

25.  Effective  slope  of  the  reef.  The  remaining  secondary  stability 
factor  is  a  combination  of  the  first  two.  This  factor,  referred  to  as  the 
effective  slope  of  the  reef,  is  obtained  by  dividing  the  cross-sectional  area 
by  the  square  of  the  crest  height.  Two  effective  slope  variables  will  be  dis¬ 
cussed  in  this  report:  (a)  the  effective  slope  of  the  structure  "as  built," 
defined  as 

A 

C’  =  ~ ^  (4) 

h'^ 

c 

and  (b)  the  response  slope  for  the  reef  breakwater  to  wave  action,  defined  as 

A 

C  =  -|  (5) 

hz 

c 

These  variables  are  considered  a  cotangent  function  since  dividing  A^  by  h^ 
one  time  produces  a  variable  which  can  be  regarded  as  a  horizontal  length,  and 
dividing  this  length  by  h^  creates  a  cotangent-like  variable.  For  low- 
crested,  or  submerged  reefs,  these  variables  provide  a  simple  way  to  charac¬ 
terize  an  average  slope  or  shape  for  what  is  sometimes  a  rather  complex  shape 
(e.g.,  see  Figure  3).  Table  4  shows  that  the  average  values  of  the  effective 
structure  slope  "as  built"  are  in  a  relatively  narrow  range.  Since  the  land¬ 
ward  and  seaward  faces  of  the  reef  were  built  to  a  slope  of  IV  on  1.5H  (cot  0 
=  1.5),  the  difference  between  the  values  of  C'  in  Table  4  and  1.5  result 
from  the  crest  width  of  the  trapezoid  which  increases  the  effective  slope,  as 
illustrated  in  Equation  6.  The  "as  built"  cross  section  of  the  reef  is  a  nar¬ 
row  trapezoid  with  a  crest  width  three  stone  diameters  wide.  For  this  study 


the  cross-sectional  area  of  the  reef  is  given  approximately  by 


■  (>■;) 
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where  cot  0  is  the  cotangent  of  the  angle  0  between  the  "as  built"  sea¬ 
ward  and  landward  breakwater  slopes  and  the  horizontal.  If  the  severity  of 
wave  attack  exceeds  a  value  of  the  spectral  stability  number  of  about  six,  the 
reef  deforms.  A  convenient  method  to  quantify  the  deformation  is  to  use  effec¬ 
tive  response  slope  for  reef  breakwaters  defined  by  Equation  5.  In  Figure  10 
the  response  slope  C  is  plotted  as  a  function  of  N*  .  This  figure  is  simi¬ 
lar  to  Figure  14.17  presented  by  Wiegel  (1964)  showing  the  relationships  among 
the  grain  size,  beach  slope,  and  severity  of  the  exposure  of  a  beach  to  wave 
action. 

26.  Because  of  the  narrow  range  of  the  effective  "as  built"  reef  slope 
C'  (Table  4),  it  was  not  possible  to  quantify  the  influence  of  this  variable 
on  stability.  It  is  assumed  that  the  flatter  the  initial  slope  of  the  reef 
the  more  stable  it  will  be.  Future  laboratory  tests  may  expand  the  range  of 
this  variable  so  that  the  influence  of  the  initial  slope  can  be  determined 
definitively. 

27.  Figure  10  suggests  that  a  logical  form  for  a  reel  breakwater  sta¬ 
bility  equation  would  be 


=  exp 


(v;) 


where  is  a  dimensionless  coefficient.  Regression  analysis  was  used  to 

determine  the  value  of  C,  for  tests  where  N*  >  6.0  ;  the  value  obtained  was 

I  s 

Cj  =  0.0945.  With  this  value  of  ,  Equation  7  explains  about  99  percent  of 

the  variance  in  C  for  the  109  stability  tests  with  N*  >  6.0  .  Equation  7 

s 

approaches  logical  limits  with 


C  ,  as  N*  -*■  °° 
s 

C  -  1 .0  ,  as  N*  -*■  0 
s 


m 

wv 


a  /ci 
Q  ?  a 


eI  ✓!' 

iH 


f  OU/5  r/OW  7 


APPROX  THRESHOLD 
OF  STONE  |  I 

MOVEMENT  \ -a  , 


Lefl  D  9 

r  □  i 


|  □  OBSERVED  DATA 


'4  16 


SPECTRAL  STABILITt  VjWBEW  Vi- 


Figure  10.  Reef  breakwater  response  slope  versus  the  spectral 
stability  number  for  stability  subsets  1,  3,  5,  7,  and  9 

since  the  natural  angle  of  repose  for  gravel  is  about  45  deg,  giving  C  =  1.0 
for  a  triangular  reef  cross  section  with  side  slopes  of  IV  on  1H.  Equation  7 
can  be  compared  to  the  observed  data  in  Figure  10.  It  is  surprising  that  the 
response  slope  of  the  reef,  stone  size  and  density,  and  severity  of  wave 
attack  can  all  be  linked  with  a  relation  as  simple  as  that  in  Equation  7.  It 
is  difficult  to  add  secondary  stability  variables  to  an  equation  like  Equa¬ 
tion  7  and  improve  the  ability  to  predict  the  response  slope  over  Equation  7 
very  much.  At  the  same  time  it  is  clear  from  Figure  9  that  secondary  stabil¬ 
ity  factors  have  some  influence  on  reef  stability.  After  trial  and  error  the 
following  equation  was  developed  which  includes  one  secondary  stability  vari¬ 
able  and  does  a  better  job  of  predicting  the  response  slope  of  the  reef: 


A  \  /h’ 

~  =  exp  <  N*  0.0676  +  0.0222  ~ 


where  the  relative  "as  built"  crest  height  of  the  reef  h'/d  was  added  to  an 

c  s 

equation  like  Equation  7  to  improve  the  predictive  ability.  Equation  8 
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explains  99.5  percent  of  the  variance  in  C  for  the  109  tests  with  N*  >  6  . 

28.  It  was  found  that  when  using  Equation  8  to  predict  the  relative 
crest  height  hc/dg  for  values  of  N*  near  or  below  six,  illogtcally  high 
values  could  result.  Higher  values  are  to  be  expected  since  Equation  6  was 
developed  for  tests  where  N*  >  6  and  there  was  enough  rock  movement  to  form 
an  equilibrium  reef  profile  and  not  for  wave  conditions  where  the  "as  built" 
reef  slope  was  too  stable  to  be  deformed.  Since  it  would  be  useful  to  have  a 
stability  model  which  predicts  reasonable  response  crest  heights  over  the 
entire  range  of  test  conditions,  another  stability  equation  was  developed  to 
predict  crest  heights  for  values  of  N*  <  10  .  This  range  provides  a  con¬ 
venient  overlap  with  the  range  of  Equation  8  and  allows  an  equation  to  be 
developed  which  will  be  simple  enough  to  serve  as  a  rule-of-thumb  relation  for 
zero  to  relatively  low  damage  situations.  This  equation  is  given  by 


=  exp 


.00005  (N*)- 
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Equation  9  provides  a  simple  relation  which  follows  the  trend  of  the  data 
well,  albeit  somewhat  conservatively  in  the  range  N*  <  10  as  can  be  seen  in 
Figure  11.  The  small  levels  of  damage  predicted  by  Equation  9  for  N*  <  6 
represent  settlement  and  consolidation  of  the  reef  under  wave  action  and  not 
conspicuous  stone  movement. 

29.  Equations  8  and  9  are  used  together  to  compute  the  response  crest 

height  of  the  reef  over  a  wide  range  of  wave  severity.  This  approach  will  be 

referred  to  as  the  stability  model.  The  procedure  is  to  use  Equation  8  for  N* 

>  10  and  Equation  9  for  N*  <  6.  If  we  let  the  solution  for  h  /h'  in 

s  c  c 

Equation  9  be  denoted  (h  /h’l  and  the  solution  for  h  /h *  in  Equation  8  be 

\  c  c)  c  c 


denoted 


(h‘/h0u 


,  then  the  following  equation 


'10  -  N* 
_ s 

dO  -  6 


\/h\  / N*  -  6  \  / h 

){v)  +  (ifcrs  (*? 

V  c  f  \  c 


can  be  used  in  the  transition  region  6  <  N*  <  10  to  compute  the  response 
crest  height  h^  .  To  judge  the  effectiveness  of  this  procedure, 
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Figure  11.  Crest  height  reduction  factor  versus  spectral  stability 
number  for  stability  subsets  1,  3,  5,  6,  and  9 

Figures  12,  13,  14,  15,  and  16  were  prepared  to  compare  observed  data  for  sub 

sets  1,  3,  5,  7,  and  9,  respectively,  with  synthetic  data  trends  generated  by 

the  stability  model.  Figures  12  through  16  show  h  / d  versus  N*  with 

c  s  s 

synthetic  trends  for  each  subset  generated  using  At  and  dg  from  Table  2 

and  hVd^  from  Table  4.  Values  of  tWd^  were  generated  at  integer  values 

of  N*  for  a  range  of  N*  about  the  same  as  observed  within  each  subset, 
s  s 

Synthetic  damage  trends  comprise  the  type  of  information  that  could  be  gen¬ 
erated  by  a  user  of  the  stability  model.  In  general,  synthetic  trends  follow 
observed  data  trends  very  well.  Discrepancies  between  predicted  and  observed 
values  appear  to  occur  because  the  stability  model  does  not  include  the  bulk 
number. 
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igure  14.  Comparison  of  data  and  the  synthetic  damage  trends 
generated  by  the  stability  model  for  subset  5 
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Figure  15.  Comparison  of  data  and  the  synthetic  damage  trends 
generated  by  the  stability  model  for  subset  7 
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Figure  16.  Comparison  of  data  and  the  synthetic  damage  trends 
generated  by  the  stability  model  for  subset  9 


Wave  Transmission 


30.  For  the  tests  mentioned  above  the  wave  transmission  coefficient  K 


is  defined  as 


Kt  -r 

c 


where  is  the  zero-moment  transmitted  wave  height,  and  is  the  zero- 

moment  wave  height  at  the  transmitted  gage  locations  with  no  breakwater  in  the 
test  channel.  Although  this  is  not  the  most  commonly  used  definition  of  K^_  , 
it  has  some  advantages  over  the  traditional  definition  which  is  given  by  the 
ratio  of  transmitted  to  incident  wave  height.  Equation  11  can  be  stated  as 
the  ratio  of  the  transmitted  wave  height  to  the  wave  height  which  would  be 
observed  at  the  same  location  without  the  breakwater  in  the  channel.  This 
definition  eliminates  wave  energy  losses  occurring  between  the  incident  and 
transmitted  gages  in  the  absence  of  a  breakwater  in  the  testing  channel. 


SS 
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These  losses  were  observed  to  be  considerable  for  the  most  severe  wave  condi¬ 


tions  during  calibration  of  the  channel.  In  effect,  measures  attenuation 


of  wave  energy  because  of  the  presence  of  the  breakwater  and  eliminates  addi¬ 
tional  energy  losses  caused  by  natural  wave  breaking  processes  occurring  be¬ 
tween  the  incident  and  transmitted  wave  gages.  Using  the  above  definition  of 


K  will  allow  evaluation  of  wave  energy  dissipating  characteristics  of  reef 


breakwaters  in  the  next  section.  Because  of  the  definition  used,  should 


be  somewhat  conservative,  i.e.,  higher  than  the  more  traditional  definition  of 


the  transmission  coefficient. 


31.  Wave  transmission  has  proved  to  be  a  very  difficult  characteristic 
of  reef  breakwaters  to  predict  partly  because  this  study  includes  both  sub¬ 
merged  and  nonsubmerged  rubble-mound  structures.  Seelig  (1980)  found  that  the 


relative  freeboard  parameter  F/H  was  the  most  important  variable  in  ex- 

mo 


plaining  wave  transmission  of  submerged  and  overtopped  breakwaters,  where 


freeboard  F  is  equal  to  crest  height  minus  water  depth,  i.e.,  F  =  -  d^  . 


However,  a  confusing  trend  will  be  obtained  using  this  variable  when  there  is 
a  transition  in  the  dominant  mode  of  transmission  from  that  due  to  wave  runup 
and  overtopping  to  that  due  to  transmission  through  the  structure.  Figure  17 


identifies  the  dominant  mode  of  transmission  as  a  function  of  the  relative 


freeboard  and  shows  a  schematized  data  trend.  The  difficulty  in  parameteriz¬ 
ing  the  wave  transmission  process  can  be  appreciated  partly  by  considering  the 
influence  of  the  wave  height.  When  a  reef  breakwater  is  submerged,  the  pri¬ 
mary  mode  of  transmission  results  from  wave  propagation  over  the  crest  and. 


generally,  the  smaller  the  wave  the  greater  the  Kt  .  When  the  crest  is  just 


above  the  water  level,  the  dominant  mode  of  transmission  results  from  wave 


runup  and  overtopping,  and  the  larger  the  wave  the  larger  the  .  If  the 


relative  freeboard  is  greater  than  about  one,  the  dominant  mode  of  transmis¬ 
sion  is  through  the  structure;  and  the  smaller  the  wave  the  greater  the  . 
A  number  of  other  factors  tend  to  further  confuse  the  above  generalities. 

32.  The  easiest  way  to  discuss  development  of  a  general  wave  transmis¬ 
sion  model  for  reef  breakwaters  is  to  first  consider  relatively  high  struc¬ 


tures  where  relative  freeboard  F/H  is  greater  than  one.  When  the  dom:'!.  nt 

mo 


mode  is  wave  transmission  through  the  reef,  is  a  function  largely  of  one 


variable  which  is  the  product  of  wave  steepness  and  bulk  number.  Figure  18 

shows  a  plot  of  K  versus  the  reef  transmission  variable  (L  d"  ^  /(H  A  ) 
t  V  P  50/  /  mo  t 


for  the  37  tests  where  F/H  >  1.0  .  This  one  variable  caused  the  wave 
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Figure  17.  Conceptual  sketch  showing  the  dominant  modes  of  wave  trans¬ 
mission  for  a  reef  as  a  function  of  the  relative  freeboard 
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Figure  18.  Wave  transmission  coefficient  as  a  function  of  the  reef 
transmission  variable  to  illustrate  the  ability  of  Equation  12  to 

predict  transmission  of  relatively  high  reefs  (F/H  >  1.0) 
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transmission  data  to  coalesce  into  one  well-defined  trend.  A  prediction  equa¬ 
tion  was  fit  to  the  data  shown  in  Figure  18,  and  the  following  relation  was 
obtained : 


1.0 


1.0  + 


0.592 


(12) 


for 


F 

H 

mo 


> 


1.0 


Equation  12  explains  about  97  percent  of  the  variance  in  for  the  range 

considered.  It  is  apparent  from  the  composition  of  Equation  12  why  the  rela¬ 
tive  freeboard  F/H  was  not  a  good  variable  for  explaining  wave  transmis¬ 
sion  through  relatively  high  breakwaters. 

33.  For  conditions  where  transmission  is  not  dominated  by  wave  energy 

propagating  through  the  reef,  relative  freeboard  F/H  is  the  most  influen- 

mo 

tial  variable.  Part  of  the  value  of  the  variable  is  in  being  able  to  account 
for  the  changing  Influence  of  wave  height  as  the  dominant  mode  of  transmission 
shifts  between  wave  propagation  over  the  crest  to  wave  runup  and  overtopping. 
For  submerged  reefs  the  relative  freeboard  correctly  indicates  the  interesting 
property  of  being  able  to  dissipate  energy  of  large  waves  more  effectively 
than  that  of  small  waves.  For  reefs  being  overtopped,  the  relative  freeboard 
correctly  indicates  that  larger  waves  have  higher  transmission  coefficients. 

In  spite  of  these  assets,  wave  transmission  for  low  and  s 'bmerged  reefs  is  far 
too  complicated  to  be  formulated  adequately  in  terms  of  F/H  ^  alone  partlv 
because  wave  energy  is  still  propagating  through  low  and  submerged  reefs  even 
though  transmission  may  be  dominated  by  either  overtopping  or  propagation  over 
the  crest.  In  addition,  energy  going  over  the  reef  is  quite  dependent  on 
crest  width  and  bulk  of  the  structure  which  introduces  the  influence  of  other 
variables.  Considering  the  multitude  of  confusing  influences  and  rfie  complex¬ 
ity  of  the  phenomenon,  the  following  regression  relation  was  fit  to  the 
167  tests  with  relative  freeboards  less  than  one: 

31 


where 


C  =  1.188 
C2  -  0.261 
C3  »  0.529 
C,  -  0.00551 

4 

Equation  13  explains  about  92  percent  of  the  variance  in  for  the  167 

tests  where  F/K  <  1.0  .  Equation  13  is  the  result  of  a  considerable  amount 
mo 

of  trial  and  error  effort  to  find  an  equation  which  fits  the  data  well,  makes 
physical  sense  based  on  current  understanding  of  the  phenomenon,  approaches 
the  correct  limiting  values,  and  is  reasonably  simple.  The  regression  analy¬ 
sis  for  Equation  13  is  shown  in  Appendix  B. 

34.  If  Equations  12  and  13  are  used,  the  transmission  coefficient  can 

be  predicted  over  the  entire  range  of  conditions  tested  in  this  study.  Pre¬ 
dicted  values  of  were  made  using  Equation  12  for  F/H  >  1.0  and  Equa¬ 

tion  13  for  F/H  <  1.0  .  This  prediction  method  will  be  referred  to  as  the 

mo 

wave  transmission  model.  Figures  19,  20,  21,  22,  and  23  show  predicted  and 

observed  values  of  K  as  a  function  of  F/H  for  subsets  1  and  2,  3  and  4, 

t  mo 

5  and  6,  7  and  8,  and  9  and  10,  respectively.  Figures  19  through  23  indicate 
that  the  wave  transmission  model  does  a  good  job  of  predicting  individual  test 
results  and  produces  trends  very  similar  to  those  of  the  observed  data. 

35.  In  addition  to  investigating  the  attenuation  of  wave  energy  passing 
ver  and  through  the  reef,  it  is  also  possible  to  determine  the  relative  shift 

:n  wave  energy  caused  bv  the  structure.  The  shift  in  wave  energy  is  measured 
the  ratio  of  the  period  of  peak  energy  density  of  the  transmitted  wave  to 
the  period  of  peak  energy  density  of  the  incident  wave.  Figure  24  shows  the 
•I .lit  in  peak  period  as  a  function  of  relative  freeboard.  What  is  surprising 


RSJS 


about  this  analysis  is  that  the  reef  does  not  produce  much  shift  in  the  peak 
period  of  the  spectrum.  In  fact,  in  only  a  few  tests  was  the  shift  as  much  as 
10  percent. 


Wave  Reflection  and  Energy  Dissipation 


36.  The  method  developed  by  Goda  and  Suzuki  (1976)  to  resolve  the  wave 
spectrum  into  incident  and  reflected  components  is  the  method  used  in  this 
study  to  calculate  the  reflection  coefficient.  According  to  Goda  and  Suzuki, 
the  reflection  coefficient  is  defined  as 


where  E^_  and  E^  are  the  reflected  and  incident  wave  energy  of  the  spec¬ 
trum,  respectively. 

37.  One  variable,  the  reef  reflection  parameter,  was  found  to  be  con¬ 
spicuously  better  than  others  for  predicting  wave  reflection  and  is  formulated 
as 


This  parameter  can  be  thought  of  as  approximately  the  ratio  of  wave  length  to 

horizontal  distance  between  the  toe  of  the  reef  and  the  swl  on  the  reef. 

Since,  for  many  tests,  the  reefs  are  deformed  and/or  submerged,  the  quantity 

(A  /h  Jd  is  sometimes  only  indicative  of  this  horizontal  distance.  When  K 
t  c/  s  J  r 

is  plotted  versus  the  reef  reflection  parameter,  a  very  strong  data  trend  re¬ 
sults  (Figure  25) .  Such  a  strong  trend  seems  surprising  considering  the  com¬ 
plex  nature  of  irregular  wave  reflection  and  the  wide  range  of  conditions 
represented  in  Figure  25.  A  regression  equation  was  fit  to  the  data  shown  in 
Figure  25  to  provide  a  convenient  rule-of-thumb  method  to  estimate  reflection 
from  a  reef  and  to  provide  insight  relating  to  wave  reflection  from  coastal 
structures  in  general.  The  equation  is  given  by 


recent  years  has  it  been  possible  to  quantify  this  property.  Quantification 
of  energy  dissipation  by  a  reef  is  the  property  that  justified  consideration 
of  rubble-mound  construction  since  both  wave  reflection  and  transmission  are 
usually  undesirable.  The  basic  conservation  of  energy  relation  for  rubble 
structures  can  be  written  as  follows: 

Kt  +  Kr  +  dissiPation  =  1*0  (15) 

where  dissipation  in  Equation  15  refers  to  the  fraction  of  the  incident  wave 
energy  dissipated  by  the  structure. 

39.  The  following  regression  equation  will  provide  an  accurate  estimate 
of  wave  reflection  from  a  reef  breakwater: 


K  =  exp 


AH 


mo ; 


(16) 


where 


C2  =  -6.77A 


-0.293 


C3  =  -0.0860 


+0.0833 


Equation  16  explains  about  99  percent  of  the  variance  in  for  the  20A 


tests  considered.  The  dependent  variables  and  the  signs  of  their  coefficients 
are  consistent  with  current  understanding  of  wave  reflection.  All  the  depen¬ 
dent  variables  in  Equation  16  affect  reflection  in  a  monotonic  manner  such 
that,  other  factors  being  equal,  K 


increases  as  d  /L  decreases,  h  /d 
2  ..  s  p  c  s 

increases,  A  /h  decreases,  and  F/H  increases.  However,  some  care 
t  c  mo 


should  be  exercised  in  using  Equation  16;  for  example,  reflection  will  in¬ 
crease  with  increasing  crest  height  only  until  the  crest  height  approaches  the 


limit  of  wave  runup  which  for  a  reef  would  be  F/H  >  1.5  .  Since  all  terms 
r  mo  ~ 


in  Equation  16  are  negative  for  submerged  reefs,  the  equation  approaches  the 


correct  limiting  value  of  =  0  for  decreasing  structure  height.  On  the 


other  hand.  Equation  16  was  fit  to  a  data  set  where  reflection  was  strongly 
correlated  to  height  of  the  reef  which  suggests  that  the  equation  might  not  be 
satisfactory  for  reefs  with  crest  heights  above  the  limit  of  runup.  This 


problem  Is  demonstrated  In  Figure  26a  where  the  difference  between  predicted 

K  and  observed  K  are  plotted  versus  relative  freeboard  F/H  .  Fig- 
r  r  mo 

ure  26a  shows  that  Equation  16  predicts  usually  to  within  ±0.05  with  lit¬ 

tle  systematic  error  except  for  high  relative  freeboards,  i.e.,  F/H  >  2.5  . 

mo 

Because  of  the  possibility  of  systematic  error  for  high  relative  freeboards, 
it  is  recommended  that  if  the  relative  freeboard  exceeds  2.5,  a  value  of  2.5 
be  used  in  Equation  16.  When  this  procedure  is  applied  to  the  data  of  this 
study,  it  removes  the  systematic  error  as  shown  in  Figure  26b. 

40.  It  is  intended  that  the  prediction  equation  for  ,  Equation  16, 
be  used  with  the  wave  transmission  model  (discussed  in  paragraph  34)  in  the 
energy  conservation  relation  given  by  Equation  15  to  compute  energy  dissipated 
by  the  reef.  This  approach  was  used  to  prepare  Figure  27  which  shows  a  scat¬ 
ter  plot  of  predicted  energy  dissipation  versus  "observed”  energy  dissipation 
caused  by  the  reef.  Figure  27  shows  that  the  procedure  outlined  above  can 
make  good  predictions  of  energy  dissipation  and  the  rather  surprising  fact 
that,  for  some  conditions,  the  reef  can  dissipate  up  to  90  percent  of  incident 
wave  energy.  Generally,  greatest  energy  dissipation  was  observed  for  short- 
period  waves  on  reefs  which  were  high  enough  not  to  be  overtopped.  The  lowest 
observed  energy  dissipation  of  about  30  percent  occurred  for  the  few  reefs 

with  a  relative  crest  height  less  than  0.7,  i.e.,  h  /d  <  0.7  .  For  sub- 

c  s 

merged  reefs,  energy  dissipation  increases  with  increasing  steepness  H  /L 

mo  p 

and  with  increasing  relative  reef  width  A  /d  L  .  Reefs  with  their  crest 

t  s  p 

near  the  swl  will  dissipate  between  about  35  to  70  percent  of  incident  wave 
energy,  and  dissipation  is  strongly  dependent  on  relative  reef  width  as  shown 
in  Figure  28.  For  reefs  with  moderate  to  heavy  overtopping,  i.e., 

0  <  F/Hmo  <  1.0  ,  energy  dissipation  is  strongly  dependent  on  the  relative 
reef  width  but  not  on  wave  steepness. 

41.  Since  wave  energy  dissipation  characteristics  of  reef  breakwaters 

are  so  important,  a  special  analysis  was  conducted  to  illustrate  the  influence 
of  the  most  important  variables  in  a  simple  way  that  would  still  be  consistent 
with  the  data.  This  analysis  used  the  most  effective  two  variables  in  pre¬ 
dicting  Kt  and  the  two  most  effective  variables  for  predicting  with  the 

provision  that  one  of  the  variables  be  common  to  both  K  and  K  so  that 

t  r 

the  predicted  values  could  be  plotted  on  a  common  axis.  Fortunately,  the  re¬ 
lative  crest  height  hc/dg  provides  a  good  common  variable.  Good  predictions 

are  obtained  for  transmission  using  the  variables  h  /d  and  B  and  for 

c  s  n 


39 


ERROR, 


wave  reflection  when  hc/dg  an<*  relative  depth  dg/Lp  are  used.  Regression 
analysis  was  used  to  develop  the  curves  for  and  shown  in  Figure  29. 

The  equations  used  to  compute  the  curves  in  Figure  26  explain  about  82  percent 
and  98  percent  of  the  variance  in  Kt  and  ,  respectively.  Appendix  B 
gives  the  equations  used  in  Figure  29  and  other  information  related  to  the 
regression  analysis.  The  curves  shown  in  Figure  29  fit  the  general  trends  of 
the  data  quite  well.  However,  the  real  value  of  Figure  29  is  that  it  is  a 
compilation  of  information  about  wave  transmission,  wave  reflection,  and  wave 
energy  dissipation  of  reef  breakwaters.  Figure  29  is  an  improvement  over  Fig¬ 
ure  8  in  Ahrens  (1984)  because  Figure  29  is  based  on  an  analytic  model; 
whereas  Figure  8  is  based  on  subjective  curve  fitting  to  the  observed  data. 
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Figure  29.  Distribution  of  wave  energy  in  the  vicinity 
of  a  reef  breakwater 
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PART  IV:  CONCLUSIONS 


42.  This  report  summarizes  the  results  from  205  laboratory  tests  of 
reef  breakwaters  conducted  using  irregular  waves.  Findings  from  this  study 
can  be  categorized  as  follows:  (a)  the  stability  of  the  structure  to  wave 
attack,  (b)  wave  transmission  over  and  through  the  structure,  (c)  wave  reflec¬ 
tion  from  the  structure,  and  (d)  energy  dissipation  by  the  structure.  These 
findings  are  largely  summarized  through  the  use  of  equations  fit  to  the  data 
which  can  be  used  to  predict  various  breakwater  characteristics  with  sur¬ 
prisingly  high  accuracy. 

43.  The  important  conclusions  from  this  study  are: 

a.  A  stability  number  was  defined  by  Equation  2  and  named  the 
spectral  stability  number  which  was  found  to  be  the  single  most 
important  variable  influencing  the  stability  of  reef 
breakwaters. 

b.  There  is  very  little  stone  movement  or  damage  for  spectral  sta¬ 
bility  numbers  less  than  six,  but  stone  movement  and  damage  can 
be  clearly  seen  for  values  greater  than  eight. 

c.  For  values  of  the  spectral  stability  number  above  six,  the  in¬ 
fluence  of  other  variables  on  stability  can  be  identified. 

Other  factors  being  equal,  the  stability  of  the  reef  increases 

the  lower  the  relative  crest  height  h  /d  ;  as  its  size  de- 

c  s 

fined  by  Equation  3  increases;  and  as  the  slope  of  the  struc¬ 
ture,  as  defined  either  by  Equation  4  or  5,  gets  flatter. 

d.  Wave  transmission  over  and  through  a  reef  is  a  very  complex 
process.  Part  of  the  complexity  relates  to  the  confusing  in¬ 
fluence  of  some  variables;  e.g.,  for  breakwaters  with  positive 
freeboards  transmission  over  the  reef  is  directly  proportional 
to  wave  height,  while  energy  transmitting  through  the  reef  is 
inversely  proportional  to  the  wave  height.  For  conditions 
where  transmission  is  dominated  by  wave  energy  propagating 
through  the  reef,  a  simple  relation,  Equation  12,  was  found  to 
predict  the  transmission  coefficient  very  well.  When  the  domi¬ 
nant  modes  of  transmission  resulted  from  wave  overtopping  or 
wave  propagation  over  the  crest  of  a  submerged  reef,  a  rather 
complex  relation,  Equation  13,  was  required  to  make  reasonable 
estimates  of  transmission  coefficients. 

e.  Wave  reflection  is  easier  to  predict  than  either  stability  or 
wave  transmission.  A  simple  relation  using  only  one  variable. 
Equation  14,  was  able  to  explain  about  80  percent  of  the  vari¬ 
ance  in  the  reflection  coefficients.  A  more  complex  relation. 
Equation  16,  was  developed  which  explained  about  99  percent  of 
the  variance  in  the  reflection  coefficient.  Other  factors 
being  equal,  reflection  coefficients  increase  with  increasing 
wave  length  and  increasingly  steeper  reef  slopes.  Reflection 


coefficients  also  increase  with  increasing  relative  reef  height 

h  /d  and  increasing  relative  freeboard  F/H  until  the 
c  s  mo 

crest  height  reaches  the  upper  limit  of  wave  runup. 


Wave  energy  dissipation  characteristics  of  a  reef  are  difficult 
to  summarize  briefly  because  of  the  complexity  of  the  phenome¬ 
non.  One  surprising  finding  was  that  for  short-period  waves 

d  /L  >  0.12  which  do  not  overtop  the  crest  the  reef  will  dis- 
s  p 

For  reefs 


sipate  80  to  90  percent  of  incident  wave  energy 

with  the  lowest  relative  crest  height  tested  0.63  <  h  /d 

c  s 


<  0.70  ,  the  structure  would  dissipate  about  30  percent  of 
incident  wave  energy.  Reefs  with  their  crests  near  the  still- 
water  level  will  dissipate  between  30  to  70  percent  of  incident 


wave  energy  depending  on  the  relative  reef  width  A  /d  L 

t  s  p 


The  model  developed  in  this  study  was  found  to  make  good  esti¬ 
mates  of  energy  dissipation. 
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APPENDIX  A:  TABULAR  SUMMARY  OF  STABILITY  AND  PERFORMANCE  DATA 
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APPENDIX  B:  REGRESSION  ANALYSIS  USED  TO  DEVELOP  FIGURE  29 
SHOWING  ENERGY  DISTRIBUTION  IN  VICINITY  OF  REEF 
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1.  For  the  energy  dissipation  figure  (Figure  29)  the  following  equation 
was  used  to  predict  the  wave  transmission  coefficient: 


•°  +  c.(r)  (Bn)C 


where 


Cx  =  0.02945 
C2  =  3.329 
CL  =  0.585 

R2  -  0.859 
F  =  611 

2.  The  wave  reflection  curves  shown  in  Figure  29  were  calculated  using 


the  following  equation: 


Kr-ex<Mci  r  + 


where 


Cj  -  0.2899 
C2  =  -0.7628 
C3  -  -7.3125 

R2  -  0.984 
F  -  4,175 
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Cross-sectional  area  of  breakwater  (cm  ) 

Bulk  number,  defined  by  Equation  3 

Response  slope  of  reef  to  wave  action,  defined  by  Equation  5 

Effective  "as  built"  reef  slope,  defined  by  Equation  4 

Dimensionless  coefficient 

Water  depth  at  toe  of  breakwater  (cm) 

1  /3 

(W-q/w^)  ,  typical  dimension  of  the  median  stone  (cm) 

h£  -  dg  ,  freeboard  of  structure  which  for  reef  can  be  either  positive 
or  negative  (cm) 

Crest  height  of  breakwater  after  wave  attack  (cm) 

Crest  height  of  breakwater  "as  built"  (cm) 

Zero-moment  wave  height  at  transmitted  gage  locations  with  no 
breakwater  in  channel  (cm) 

Zero-moment  transmitted  wave  height  (cm) 

Incident  zero-moment  wave  height  (cm) 

Reflection  coefficient  of  breakwater  as  defined  and  calculated  by 
method  of  Goda  and  Suzuki  (1976) 

H  /H  ,  wave  transmission  coefficient 
t  c 

Airy  wave  length  calculated  using 

Stability  number,  defined  by  Equation 

Spectral  stability  number,  defined  by 

Wave  period  of  peak  energy  density  of 

3 

Density  of  stone  (g/cm  ) 

3 

Density  of  water,  tests  conducted  in  fresh  water,  ww  ”  1-0  (g/cm  ) 

Median  stone  weight  (subscript  indicates  percent  of  total  weight  of 
gradation  contributed  by  stones  of  lesser  weight)  (g) 


and  d  (cm) 
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